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Equation of state parameter plays a significant role for guessing the real nature of dark energy. 
In the present paper polytropic equation of state p — ujp" is chosen for some of the kinematical 
A-models viz., A ~ (a/a)'^, A ~ a/a and A ~ p. Although in dust cases {uj — 0) closed form solutions 
show no dependency on the polytropic index n, but in non-dust situations some new possibilities 
are opened up including phantom energy with supernegative (a; < —I) equation of state parameter. 

PACS numbers; 04.20.-q, 04.20.Jb, 98.80.Jk 



I. INTRODUCTION 

Ever since the discovery of an accelerating Universe 
through SN la observations P, 01 , scientists are search- 
ing for the cause behind this acceleration. It has been 
suspected that some kind of yet unknown energy is re- 
sponsible for injecting the right amount of energy to the 
Universe for changing it to an accelerating one from a 
decelerating phase. Scientific community has coined the 
name dark energy for this unknown energy. Various types 
of models have been proposed for approximating this 
dark energy so far 0, Ql • 

It is observed that equation of state parameter plays 
a crucial role in understanding the actual nature of dark 
energy Till now, most of the models of A (the so- 
called cosmological constant of Einstein and a represen- 
tative symbol for dark energy) with dynamic character 
have relied on the barotropic equation of state with wide 
range of values of the equation of state parameter uj, viz., 
u! = (for dust filled Universe) , ui — 1/3 (for radiation) , 

= — 1 (for vacuum), — 1 < w < (for quintessence) 
and UJ < ~1 (for phantom energy). But, all these models 
are plagued by some shortcomings. For instance, models 
with Lu — and 1/3 show an excess 6] or very low 3] 
age of the Universe while the stiff-fluid model, in spite of 
its nice agreement with the present age of the Universe 
0, is not, in general, accepted as the real nature of the 
present Universe. In this circumstances, it is not unrea- 
sonable to think of an equation of state different from the 
barotropic one. This has prompted us to investigate the 
reaction of some of A models when polytropic equation 
of state is chosen. 

Polytropic equation of state has been used in various 
astrophysical situations such as in the case of Lane- 
Emden models H, 1^ . In the present investigation it is 
used in cosmological realm. As test models we have 
selected A ~ [a /of', A ^ d/a and A ^ p models - the 
same three kinematical A-models which were chosen 
in one of our previous works |lClj | under the barotropic 



equation of state. Here Sections 2 and 3 deal with the 
fleld equations and their solutions under three different 
A-models while various physical implications of the 
present work are discussed in Section 4. 



II. EINSTEIN FIELD EQUATIONS 



The Einstein field equations are 



2 ^ 



-SttG 



A 
8^ 



9- 



(1) 



where cosmological constant A is assumed as a function of 
time, viz., A — A{t) and the velocity of light c in vacuum 
is unity when expressed in relativistic units. 

For the spherically symmetric Friedmann-Lemaitre- 
Robertson- Walker (FLRW) metric 



ds^ = -dt^ + a{tf 



dr^ 



1 — fcr^ 



+ r^{d0'^ + . 
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where a is the scale factor and fc, the curvature constant 
— 1,0,-1-1 respectively for open, flat and close models of 
the Universe, the Einstein field equations (1) take the 
forms as 
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As various observational results [ll| and inflation theory 
indicate that the Universe is flat so we have assumed 
here k = Q. 

Let us choose the polytropic equation of state 

p = wp" (5) 
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where equation of state parameter w can take the 
constant values 0, 1/3, —1 and +1 respectively for the 
dust, radiation, vacuum fluid and stiff fluid and n is the 
polytropic index. 



III. PARTICULAR SOLUTIONS 

A. A ~ {ii/af 

If we use the ansatz 

A = 3«g)' (6) 

where a is a constant, then using equation (6) we get 
from equation (3) 



3 ( ^ y _ 3a y = SttGp 



which takes the form for density as follows: 



3(1 ~ a) f a 
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Using equation (5) and equation (6) we get from equation 
(4) 



3 ( ^ j = -127rGcjp" -47rGp + 3a f ^ ) (9) 



which, on simplification, yields 



_i (Pa ^ f da 



2n 



da 



where A = [3"w(l - a)"]/[2"(47rG)"-i] and B = [1 - 
3a)/2. 

Let us substitute, z — da/dt for which the equation 
(10) reduces to 



da a 



(11) 



If we put, ^" = J/, then differentiating it with respect 
to the scale factor a, we get 



i-2ndz _ 1 dy 



da (2 — 2n) da 
Hence equation (11), by use of equation (12) becomes, 
dy (2 - 2n)B 



(12) 



da a 



-y -(2 - 2n)Aa^ 



(13) 



Now let us solve equation (13) in some particular cases 
for getting insight into the physical situations. 



1. Dust case (u — 0) 

For pressureless dust, uj — — A and equation (13) 
becomes, 



dy (2 - 2n)B 



da a 
Solving equation (14) we get, 

y = Gia(2"-2)S 

Gi being an integration constant. 
This immediately yields 

da _^ B 



y = o. 



(14) 



(15) 
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where G2 = Gi 1/(2-2"). 

Integrating equation (16) and using the initial condition 
that a = for < = we get. 



i{t) = Gs 



3(1 - a) 



2/3(1-q) 



^2/3(l-a) 
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where G3 — C2^^''^^^^ being an integration constant, 
which provides the Hubble parameter as 



H{t) 



a 
a 



3(1 -a)t' 



Therefore from (8), we get 
Pit) 



1 



1 



67rG(l - a) t^' 
Also, from equation (6), we obtain 

4a 1 



3(l-a)2i2- 



(18) 



(19) 



(20) 



Thus, it is interesting to note that in dust case, ex- 
pressions for p{t) and A{t) as well as t dependent part 
of a{t) are independent of n. Moreover, a{t), pit) and 
K{t) follow the same power law with t as those in the 
dust case with barotropic equation of state (i.e. n = 1) 
10]. But, surprisingly, equation (17) becomes undefined 
for n = 1 since G3 contains the factor (2 — 2n) in the 
denominator. 



2. Non-dust case (uj ^ Q) 

When a; ^ 0, then A^Q and hence equation (11) is a 
linear equation. Multiplying equation (11) by integrating 
factor a(2-2n)-B solving the resulting equation we get. 



Aa^ 



2-2n)B 



B+l 



(21) 
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which yields after substituting the value of y 
da A 



dt (B + 1) 



(22) 



By solving equation (22) we get our solution for the scale 
factor in the form 

a{t) = C4e.p[(-c.)V2(i-)^^!I^]i/2i (23) 

3{1 — a) 

where C4 is integration constant. 

It is clear from equation (23) that no real value of a{t) is 
possible if to is positive. So, uj must be negative. Putting 
specifically, w = — 1 in equation (23) we obtain 



a{t) = CiexpiSnG/iil - a)]^^^t. 



(24) 



Now, Lu = —1 means a vacuum fluid. So, for vacuum 
fluid, we get an exponential solution independent of n. 
But, ii u! < —1, then scale factor depends on n. Now, 
cij < — 1 corresponds to a supernegative equation of state 
and hence the idea of phantom energy comes into the pic- 
ture. Also, during inflation, the Universe underwent an 
exponential expansion. So, for non-dust case our solution 
reflects an inflationary scenario which may be due to ei- 
ther vacuum energy or quintessence or phantom energy. 
Also, from equation (24) we get, 



H = 



_8ttG_ 

3(1 - a) 



1/2 



(25) 



From equation (8) and (6) we get respectively 

p{t) = 1, (26) 



m = 



SnGa 
1-a' 



(27) 



Therefore, the cosmic matter and vacuum energy densi- 
ties are given by 



nr. 
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So that, equations (29) and (30), immediately provide 



rim + f^A = 1- 



(30) 



Thus, we find that although polytropic equation of state 
yields interesting result by invoking the idea of phantom 
energy for non-dust case, yet it presents us some unac- 
ceptable situations like constant H (equation (25)) and 
constant p (equation (26)). A possible explanation of 
this will be discussed afterwards. However, the cosmic 
matter and vacuum energy density parameters satisfy 
the well-known relation Qm + — 1 for flat Universe 

{n^ = 0). 



B. A ~ p 

For this model we use the ansatz 

A = &t:Gjp 

where 7 is a constant. 

Then, from equation (3) we get 



87rG(7 + l) \a 
Using (5) we get from (4) 
'a 



(31) 



(32) 



= -U-kGljp'' - AiiGp + 87rG7p (33) 



which, on simpliflcation, yields the differential equation 

^ + f ^ V" _ f = (34) 

dt^ \dt) \dt) ^ ' 

where D = w3"/[2"(7 -I- l)"(47rG)"-i] and S = (27 - 
l)/2(7+l). 

Substituting da/dt — u and then v — u'^ ^" equation 
(34) becomes 



1 dv E 



(2 — 2n) da a 



V = -Da' 



(35) 



1. Dust case (uj — 0) 

For pressureless dust case D—0, so that equation (35) 
reduces to 



dv , ^ E 

— + (2n-2)-v = 0. 

da a 



(36) 



Solving equation (36) in the same manner as in the pre- 
vious model, we get our solution set as 



ait) = C5 



2(7 + 1) 



2(7+l)/3 



t 



where C5 is a constant. 



Hit) 



Pit) = 



2(7+1)1 



(7 + 1) 1 

67rG t2' 

47(7+1) 1 
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(37) 



(38) 



(39) 



(40) 



Also, we can calculate for the cosmic and vacuum energy 
densities which, respectively, are Qm = 1/(7 + 1) find 
^A = 7/(7+1)- Therefore, + 51 a — 1 and hence 
7 — ilA/^m- Thus, we find that for this model also ait), 
pit) and A(i) follow the same relationship with time as 
in the barotropic case 10] and 7 is related to matter and 
vacuum energy densities by the same relation as in the 
barotropic case p^ . 
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2. Non-dust case (uj ^ Q) 

When (jj ^ 0, then D 0. In this case, following the 
same procedure as in the A ^ {a /of' model, it is easy to 
show that no real solution is possible for w > 0. 
For a; = — 1, we have 



a{t) = Ceexp 
where Cg is a constant. 



B7rG(7 + l) 



1/2 



(41) 



Here we find that in non-dust case, solutions are possible 
for vacuum fluid {lo = —1), quintessence (—1 < uj < 
0) and phantom energy (w < —1). For vacuum fluid, 
solution is independent of n whereas for quintessence and 
phantom energy solution depends on n. In this case also, 
7 is found to be the ratio of vacuum energy density and 
matter energy density. For physical reality, 7 > — 1 to be 
imposed on the solutions. 



C. A ~ a/a 



H 



Pit) 



3 



-,1/2 



87rG(7 + 1) 
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87rG(7 + 1) 



2 " 



(42) 



(43) 



(44) 



If we use the supposition A = (i{a/a) = f3{H + i?^), 
then from equation (3) we get, iH'^ 87rG/9+/3(ij+i7^). 
Therefore 



1 



87rG 



;(3-/3)ij2 -/3ij]. 



(45) 



Then, using (5), from (4) we have, Z{H + H"^) = 
— 127rGwp" — 47rG/9 + (3{H + iJ^) which on simplification 
becomes 



(3 - /3)if2 + (2 - p)H 



127rGw 
(87rG)" 



[(3 - I3)H'^ - PHY". 



(46) 



For simplifying equation (46), let us investigate in some 
particular cases. If we specifically choose n = (i.e. 
p ~ uj), then equation (46) reduces to 



(2 - I3)H + (3 - (3)H'^ = -87rGw. 



(47) 



1. Dust case (lo — Q) 
For pressureless dust [lu — 0) equation (47) becomes 

(48) 



{2-(3)— + {3-(3)H' = 0. 



By solving equation (48) and using the initial condition 
that a = when t = we get 



da 



2- (3\ 1 



(49) 



dt \3-f3j t' 
Then from equation (49) we finally get our solution set 



where G7 is a constant. 



Pit) 



1 (P-A 1 

47rG \(3-3j f^' 



(50) 



(51) 



m = 



Hit) 



m - 2) 1 

(/3-3)2 i2' 



/3-2\ 1 



P-3 J t 



(52) 



(53) 



Thus, for n = 0, the solution set is same for a(t ), p (t) 
and k{t) as obtained by Ray and Mukhopadhyay [i3 
the dust case with barotropic equation of state. 



2. Non-dust case (uj ^ Q) 



In non-dust case, equation (47) can be written as 

J TIT 

(2 -/?)— + (3 - I3)H^ = -87rGw. 
dt 



(54) 



By solving equation (47) and remembering that H — 
when t — we get 



da 
It 



va tan{pt) 



(55) 



where ^= [87rGw(3-/?)/(/3-2)]i/2 and = [87rGw/(3^ 
f3)]'/'- 
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After solving equation (55) we have 
{87rGw(3 /?)} 



1/2 T {87rGa./(3-/3)}i/^ 



(/3-2) 



,(56) 



P(i) = 



(87rGa;)3/2 (3-/3) 
87rG(/3 - 2)2 



{87rGw(3-/3)}i/2W 



,{87rGu;(3-/3)}i/2 ^ {87rGc.(3 - /3)}i/2 



(/3-2) 



— /3sec 



(/3-2) 



(57) 



A(t) 



(87rGco)3/2 
(/3-2)2 



/3 



(3-^)Vwi^!^^^4^^M^t- 



(87rGw)i/2tan' 



{8^Glj(3 - /3)}i/2 
(/3-2) ' 



(58) 



87rGu; {87rGw(3 
-tan- ^ 



/3)} 



1/2 



(/3-2) 



(/3-2) 



(59) 



From the above solution set, it is clear that if > 0, then 
for real a(i), p(t), A(i) and i?(i), /9 must be less than 
3. For negative w (i.e., for vacuum fluid, quintessence 
and phantom energy) real values of ait) and Hit) can 
be obtained, but A(t) and p{t) become imaginary. If 
(/3 — 3) = 0, then p{t), A(t) and -ff(t) are zero but a(t) 
is undefined. Since tant and sect are both increasing 
functions of i, then a(t), A(t) and Hit) increase with 
time. Increasing a(t) supports the idea of an accelerating 
Universe, but increasing A(i) and Hify are contrary to 
the present status of those two parameters. 
For i = 0, we get respectively from equations (56) - (59) 
the following physical parameters: 



p{t) = -13 



a(t) = Cs, 



1 (8^Gu;)3/2 (3-/3) 
87rG (/3 - 2)2 



(8^G^)3/2 
(/3 - 2)2 



/3(3-/3)i/2, 



H{t) = 0. 



(60) 



(61) 



(62) 



(63) 



Since a{t) assumes a constant value for t — 0, then our 
model hints at the existence of some kind of quantum 
fluctuation at the time of Big Bang. Now, equation (61) 
tells us that for physically valid p, uj > and /3 < 0. For 
a negative /3 we get an attractive A. So, as a whole this 
model also presents some interesting as well as awkward 
cosmological picture. 



IV. DISCUSSION 

Present investigation reveals that for pressureless 
dust, all the three models reflect the same result as 
obtained by Ray and Mukhopadhyay '^'^^^ '^^^e 

with barotropic equation of state. But, in non-dust 
cases (i.e., a; 7^ 0), A ~ (a/a)2 and A ~ p models show 
exponential expansion of the Universe for negative uj 
whereas no real situation is possible for lu > 0. Thus, 
we can say that for dust-filled Universe, there is no 
distinction between barotropic and polytropic equations 
of state. This is not at all unexpected because for uj = 0, 
we have p = 0, whatever may be the value of n in the 
equation of state. Also, non-dust cases for all the three 
models present us some unpleasant results in terms of 
constant H and p. These situations can be explained 
if we assume that the non-dust cases reflect the picture 
of early Universe when inflation occurred due to the 
presence of quintessence (— 1 < w < 0) or vacuum fluid 
(lu ~ —1) or phantom energy [uj < —1). So, using 
polytropic equation of state it has been possible to show 
that non-dust cases admit the presence of a driving force 
behind inflation in the form of either quintessence or 
vacuum fluid or phantom energy and in the dust cases 
there is no distinction between different equation of 
states. Moreover, present models do not depend on any 
particular value of n in the polytropic equation of state. 
This proves the generality of the present investigation. 
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